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Abstract 

We obtain the exact solution of the Klein-Gordon equation describing the 
propagation of a particle in two regions of different constant magnetic field, 
separated by an infinite plane wall. The continuity of the wave function and of 
its derivative at the interface is satisfied when including evanescent-wave terms. 
We analyze solutions on truncated spaces and compare them with previously 
obtained approximate solutions. The findings of this work have applications in 
the problem of the propagation of particles in the presence of a bubble wall in 
the midsts of an electroweak phase transition, where the two regions separated 
by the wall are influenced by different (hyper)magnetic field strengths. 
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1 Introduction 

The possible existence of magnetic fields in the early universe has recently become 
the subject of intense research due to the various interesting cosmological implica¬ 
tions mi. Among them is the effect that primordial magnetic fields could have 
had on the dynamics of the electroweak phase transition (EWPT) at temperatures 
of order T ~ 100 GeV |3J. It has been recently pointed out that, provided enough 
CP violation exists, large scale primordial magnetic fields can be responsible for a 
stronger first-order EWPT ®E! in analogy with a type I superconductor where the 
presence of an external magnetic field modifies the nature of the superconducting 
phase transition due to the Meissner effect. Primordial magnetic fields, extending 
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over horizon distances, can be generated for example in certain types of inflationary 
models through the breaking of conformal invariance |5]. 

Recall that for temperatures above the EWPT, the SU(2)/ / xU(l)y symmetry is 
restored and the propagating, non-screened vector modes that represent a magnetic 
field correspond to the U(l)y group, instead of the U(l) em group, and are therefore 
properly called hypermagnetic fields. 

Another interesting implication of the presence of these fields during a first-order 
EWPT concerns the scattering of fermions off the true vacuum bubbles nucleated 
during the phase transition. In fact, using a simplified picture of a first-order EWPT, 
it has been shown [7J that the presence of such fields also provides a mechanism, 
working in the same manner as the existence of additional CP violation within the 
SM, to produce an axial charge segregation. The asymmetry in the scattering of 
fermion axial modes is a consequence of the chiral nature of the fermion coupling to 
hypermagnetic fields in the symmetric phase. 

The simplified treatment of the scattering process allows to formulate the problem 
in terms of solving the Dirac equation with a position-dependent fermion mass jSJ. 
This dependence is due to the fact that fermion masses are proportional to the vacuum 
expectation value of the Higgs field, this last being zero in the false phase and non¬ 
vanishing in the broken symmetry phase. Transmission and reflection coefficients can 
be computed approximately by looking into the most likely transitions from a given 
incident flux of particles. 

Nevertheless, in spite of the success of such a scheme, it is clear that a more 
accurate account of such processes is needed in order to set better limits that constrain 
the feasibility of the mechanism for axial charge segregation ||. 

In this work, we present an exact solution to the problem of describing the propa¬ 
gation of a particle in two regions of different constant magnetic field, separated by an 
infinite plane wall, corresponding to the wall of a bubble nucleated during a first-order 
EWPT. The particle is massless outside the bubble, in the symmetric phase, and mas¬ 
sive inside the bubble, in the broken-symmetry phase. The analysis is made in the 
infinitely thin-wall approximation to avoid unnecessary numerical complications. In 
physical terms, this means that the spatial region over which particle masses change 
is small compared to other relevant length scales such as the particle mean free path. 
The choice of external field captures such information, namely, changes in the mag¬ 
netic field occur over large spatial regions, and thus can be considered homogeneous 
near the interphase. 

Finite temperature effects are implicitly included in the treatment of the problem. 
One could model the change of the particle’s mass with position using the so-called 
kink solution for the spatial profile of the Higgs field j8|. This profile is obtained from 
the classical solution to the equations of motion with the finite-temperature effective 
potential. The wall profile considered in this work can be thought of as obtained from 
the kink solution in the limit where the wall thickness goes to zero. 

In order to provide a working solution, we perform a numerical analysis by trun- 
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eating the solution space. To present the method, we choose to work within the 
framework of the Klein-Gordon equation although applications to the Dirac equation 
are straightforward. 

The outline of this work is as follows: In Sec. [21 we briefly review the well-known 
solution to the problem of a spinless, charged particle moving in a constant magnetic 
field. In Sec. 01 we solve the problem that describes the motion of this particle 
propagating through a planar interface dividing two regions with different magnetic 
field strengths, setting the system of equations that allows to find its reflection and 
transmission amplitudes. In Sec. 0 we find a numerical solution by truncating the 
infinite tower of coupled equations that result from the exact treatment. In Sec. 01 we 
use the above solutions to compute reflection and transmission coefficients for modes 
moving from the symmetric phase toward the broken-symmetry phase and explicitly 
check that these satisfy flux conservation. In Sec. [01 we build a one-dimensional model 
to understand such flux conservation in spite of the relativistic context of the problem. 
In Sec.|7]we summarize this work and give an outlook about its implications. Finally, 
in the appendix we explain how the proposed method is applied for arbitrary constant 
magnetic fields. 

2 Klein-Gordon equation and solutions for con¬ 
stant magnetic field 

We start by describing the motion of a particle with charge eo and mass /i, under a 
constant magnetic field of magnitude B. Working in a gauge where the time compo¬ 
nent of the corresponding vector potential is zero, we can choose 

A = _ ^ rXB ’ W 

with B = (0,0,5). By minimal substitution p — ► p + eoA in the Klein-Gordon 
equation we have 

[-V 2 + ^e 2 5 2 r 2 + e 0 BL z + /i 2 ]T = 5 2 T, (2) 

where we use h — c—1 and, working in cylindrical coordinates (r, (j),z),r= yjx 2 + y 2 
is the radial coordinate. Also, L z = —i(d y x — d x y) is the ^-component of the orbital 
angular momentum. 

The solutions are given in various representations in Ref. mini. In one choice, 
a complete set of eigenfunctions is characterized by an angular momentum quantum 
number l G Z, and we assume l > 0; a radial transverse xy -plane excitation quantum 
number s = 0,1, 2...; and free-particle propagation along the magnetic-field direction, 
with wave number k n , which with our choice of magnetic field direction as in Eq. 0 , 
points along the z axis. The energy is given by 

E = yfkl + 2 e 0 B(n + 1/2) + p 2 , (3) 
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and the eigenfunctions by 
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where iV ns = (e 0 B/E) 1/2 , p = ( e 0 B/2)r 2 , n = s + l, L is the distance along 0 and 
it eventually cancels out in the calculation of reflection and transmission coefficients. 
The definition of N ns differs from that of Ref. PI by a yU 1 / 2 factor. For k n real, the sign 
of the exponent in e lknZ represents the direction of motion of the wave. We note that 
the set of quantum numbers (n, s, k n ), (n, /, k n ), ( s , l, k n ) can be used interchangeably. 

The normalization coefficients are such that the set of eigenfunctions satisfy the 
normalization condition 


/ 
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where p 0 is the charge density. Q l s {p) is a Laguerre polynomial defined by 


QsiP ) = s\L l a {p), 


(7) 


where we use the convention 


s 


L 1 S (P ) = E 

5=0 


(s + l)\p s j 
j\{s - j)\{s + l - j)\' 


( 8 ) 


For negative /, the solutions are Q l s (p) = (—1 ) l p l Q s h(p), and thus the energy eigen¬ 
value is E = i \Jk 2 + 2 e 0 B(s + 1/2) + p?. Correspondingly, a new complete basis set 
over s on the xy plane can be defined. 

A general solution T ei can be constructed from a linear superposition of the basis 
of solutions in Eq. 0 and is given by 


^Ei{r, <f>, z) = E C n fB S (r, 0, z), (9) 

n 

where C n are arbitrary coefficients that weigh the corresponding Landau-level contri¬ 
butions. 

All such modes are in principle necessary in the context of scattering, where an 
incoming wave is either reflected or transmitted conserving energy, but the magnitude 
of its momentum eigenvalue along the 5 axis, k n , may change. Indeed, the energy 
gained or lost along the z axis is redistributed in the form of radial excitations asso¬ 
ciated with the quantum number s. Thus, changes in s are compensated by changes 
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in k n in such a way that the energy is conserved, and is equal to the one in Eq. ©• 
Notice that, for large enough n and fixed energy E, k n will eventually become imag¬ 
inary. We refer to this kind of solution as evanescent wave. Physically, this situation 
mimics the transmitted wave function found in the presence of a square barrier with 
a potential value larger than the energy E. 

That momentum along the z axis is not necessarily conserved is physically clear 
from the fact that translational invariance along the 5 axis in this problem is lost due 
to the presence of the wall, as will be exemplihed in the one-dimensional models in 
Sec. El 


3 Solution for two regions with different magnetic 
fields separated by an infinite plane wall 


We proceed to find the solution for the propagation of a particle in two regions with 
different constant magnetic fields, and different mass. We use the solutions found in 
the previous section, placing our coordinate system with the xy plane on the wall 
separating the two regions. For definiteness, we call region / the half-space with 
z < 0, where the held strength is B, the coupling constant e 0 , and y = 0, and region 
II the half-space with z > 0, where the held strength is B', the coupling constant e' 0 , 
and y ^ 0. The magnetic helds keep perpendicular to the wall. 

An incoming hux of charged, spinless and massless particles from region /, propa¬ 
gating along the z axis, scatters oh the wall at z = 0, and is rehected back into region 
/, and transmitted into region II. 

We prepare the incoming hux with each particle propagating with momentum k l m , 
quantum number s -representing a given transverse oscillation mode- and angular 
momentum l. These values define E. The particle’s wave function satishes the Klcin- 
Gordon equation with magnetic helds B and B', corresponding to regions / and II, 
respectively. Notice that the energy and angular-momentum component along the z 
direction are conserved. This means that E and l will be good quantum numbers on 
both sides of the wall. Also, we demand that the wave function and its derivative be 
continuous at the interface. 

We normalize the incoming wave function to unity, as in Eq. El which from Eq. © 
means that C^ = 1. Among the scattered waves, one corresponds to a rehected wave 
with magnitude of the momentum k m equal to the incident one, but with direction 
of motion reversed 


Ki = ~k rn = \JE 2 - 2e 0 B(m + l/2 ); (10) 

we associate to it the coefficient C^, which represents the reflection amplitude. The 
total wave function in region / can be expressed, according to Eq. © , as 
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where in the first term we have used that 
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as implied by Eq. ©• 

Recall that given an energy eigenvalue E and an orbital angular momentum l, to 
each eigenfunction-labeled by n-among the inhnite set, corresponds a value of the 
momentum k n . Thus, from Eq. ©, the possible values for k n are given by 

_ J —a/ E 2 — 2eo B(n + 1/2), E 2 > 2eoB(n + 1/2) , . 

n ~ \ —i\J—E 2 + 2 eoB(n + 1/2), E 2 < 2e 0 B(n + 1/2). 1 j 

Notice that the majority of the k n are imaginary, for there is a value n beyond which 
the square-root argument is negative. The sign of the imaginary k n is chosen so as 
to have a finite wave function as z —> — oo. 

Also, a general solution in region II can be written as 


(p,z) = J2 C n'ft'\r, 0, z), (14) 

n' 

where the solutions have the same l and E, and a quantum number n' given by n! = 
s' + l, with a value of the longitudinal momentum, chosen to represent transmitted 
or tunneled evanescent waves, given by 

r y/g - 2+ 1/2) - E 2 > 2e' 0 B'(n’ + 1/2) + jj? 

1 i^J-E 2 + 2e' 0 B'(n' + 1/2) + p 2 , E 2 < 2e'„B'(n’ + 1/2) + f, 2 , 1 ’ 

where the sign of a real k n ' represents a wave travelling in the +z direction, and the 
sign of an imaginary k n i is chosen so as to have a finite wave function for z —> oo. 

The constants C^, C„, CjJ in Eqs. (fill) and (TTH) are coefficients to be determined 
from the continuity of the solution and its derivative across the planar interphase. 

First, in order to join the solutions from regions / and II, we can re-express the 
region -II transverse wave functions in terms of the transverse-component basis of the 
region-/ solutions at z = 0: 


/S:*'(r,^,0) = ^a‘ y /g*(r>,0), 

S 

where the z dependence factors out as in Eq. ca 


/S'/(r>, Z ) = /"(r,0.O)e“»". 

The amplitude a! ss , can be calculated in a closed form as 

1 f s's' 1 \ 1 ^ 2 I°° 

<C - JAUW I jra) (e 0 e' 0 BB'/i )‘/ 2 J dxt. '"i 1 


( 16 ) 


(17) 


( 18 ) 
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where we have changed variable with x = r 2 , and have carried out the angular 
integration; we have also defined A = e 0 B/2, rj = e' 0 B'/2, b = ( e 0 B + e' 0 B')/A. The 
integral in Eq. m is given by [13] 


dxe- bx x l L l s (\x)L l s ,(r]x ) = 


T(s + s' + l + 1) (b — A ) s (b — 7]) s ' 


sis'! 


^s+s'+Z+l 


X 


(19) 


FIs' b ( b ~ X -ri ) 

1 ’ ’ ’ (b-\)(b-nr 


where T is the gamma function, and F is the hypergeometric function 2 iA(a, 6; c,z). 

Thus, in Eq. (11411 can be written in terms of the basis of in Eq. (IIID . by 
means of Eqs. m and m , as 


<t>,z) = ^2 C^a' ss ,f™(r, </>, 0)e 

n' ,s 


ik„,z 


( 20 ) 


The equality of ty T El and at z — 0 implies 


cl = Y. a '“' c "' n * 

n' 

1 + Cm — a L.s'CnM 


m 


( 21 ) 

( 22 ) 


where each n mode is matched (as l is fixed, summing over n! implies summing over 
s'.) Notice that this matching also implies that of all the transverse functions. We 
also impose the equality of the derivatives across the surface, namely, 

d^Ei/dz = d*%/dz (23) 


at z = 0, which implies 


Cik n = J2 al ss' C n' k n’, ^ + 

n' 

k m + C'rn.km = (~ 1 + C^)k m = ^ a l SmS ,CfJk n i , 


m 


(24) 

(25) 


where we also take account of the reflected wave m, with s m = m — l, that propagates 
with momentum k m , opposite to that of the incident wave. We thus find an infinite 
system of non-homogenous linear equations for CC ^ and CfJ, representing reflec¬ 
tion and transmission amplitudes, respectively. Their actual solution is analyzed in 
Section 0 This method is generalized in the appendix for magnetic fields with ar¬ 
bitrary direction. While we have assumed an incoming wave with definite quantum 
numbers, the solution is general for any incoming wave function, as one only needs 
to expand it in the basis of Eq. m, and solve for each component. Indeed, a wave 
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function |\P 0 ) representing initial conditions can be expanded in terms of the solution 
basis: 


*o> = |*o)|sifc I ), 


(26) 


Slk ; 


where an integral over the momentum k z is implicit. 

4 Truncated solutions 

In order to find the numerical solution to the system of Eqs. (EH), EH and EH, EH, 
we need to truncate it. For an incoming particle in region / with associated radial 
quantum number m, the relevant scattering states into which it can be reflected are 
those with momentum eigenvalue close to k m , and therefore with radial quantum 
numbers n close to the incident one. For the transmitted states on side //, it has 
been argued H that for similar values of e^B and e' 0 B' the main contribution comes 
from the m! = m state and its closest neighbors. Concentrating on such a case, it 
is sensible to formulate the above system of equations considering the states from 
m — n max to m + n max , and check the convergence for given integer n max . In fact, as 
the relevant m is zero or close to it, we shall consider n e [0, n max \. Such a truncation 
is appropriate since, as exemplified in Fig. |T| for C{, the coefficients quickly converge 
for a finite value of n max - 

In such a case, the matrix representation of the above system of equations can be 
given in terms of four blocks, each one of size n max x n max , 



(27) 


where the upper blocks contain Eqs. EH>, EH, and the lower ones Eqs. EH, EH 
(the notation used keeps the CC T J indexes), and the terms with <5 n » m contain the 
non-homogenous part, and define the incoming m state. 

With the method thus described, we have obtained a quick convergence of the 
reflection and transmission coefficients as function of n max , computed using a set of 
parameters corresponding to the analysis in the following section. These coefficients 
do not change after a relatively small value for n max which for the chosen set of 
parameters corresponds to n max ~ 20. The sum of the reflection and transmission 
coefficients adds up to unity. This point is further discussed in the following sections. 

5 Transmission and reflection coefficients 

To compute the reflection and transmission coefficients, we consider the a flux of in¬ 
coming positive-energy particles from region / with definite values of k l m and quantum 


c 


I 

1 



Figure 1: Sensitivity of coefficient C[ as a function of the maximum number of 
coefficients in the truncated solution, n max . The magnitude of the hypermagnetic 
held in the symmetric phase is B = 0.01 GeV 2 . The magnitude of the magnetic held 
in the broken-symmetry phase is B' = cos 9\yB. The incoming wave function has 
quantum numbers s — 0, l — 1, an incident energy E = 1.5 GeV and mass fi = 1 
GeV. The coupling constant in the symmetric phase is taken as eo = g' = 0.344, 
whereas in the broken-symmetry phase it is e’ Q = e = g'cos9w- 


numbers s and l. The conserved current is given by 

f = - eo&'&A 1 *. (28) 

We are interested in the longitudinal component of this current integrated over the 
transverse x—y plane at z —> —oo for the incident and reflected waves and at z —► +oo 
for the transmitted wave, namely 

3 = j j z dxdy, (29) 

where j represents the incident, reflected or transmitted integrated currents 

ii = *4 (30) 

j, = J]Re(fe„)|C4 2 

n 

j, = ^Re(i„,)|C';| 2 . 

n' 
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The reflection and transmission coefficients R, T , are built from the ratio of the 
reflected and transmitted integrated currents along the z direction, to the incident 
one: 


R = ~jr/ji , (31) 

T = jt/ji. (32) 

The incoming wave function corresponds to the first term on the the right-hand side 
of Eq. (1111) . whereas the reflected one corresponds to the second and third terms. The 
transmitted wave function is given by Eq. EH- The coefficients C„, C ^ and CfJ are 
found numerically from the solution to Eqs. (EH), EH and EH, EH by truncating 
the solution space as described in Sec. 

Figure [21 shows R and T for the case of an incoming wave function with quantum 
numbers s — 0, l = 1. The magnitude of the hypermagnetic field in the symmetric 
phase is B = 0.01 GeV 2 . The magnitude of the magnetic field in the broken-symmetry 
phase is obtained from the requirement that the Gibbs free energy be a minimum, 
and is given by jo 


B' = cos OwB , (33) 

where 6w is Weinberg’s angle. The coupling constant in the symmetric phase is 
associated to the hypercharge U( l)y symmetry, so that e 0 = g' = 0.344. The coupling 
constant in the broken-symmetry phase is associated to the electromagnetic U( l) em , 
that is, e' 0 — e — g' cos 9\y- The mass in the broken-symmetry phase is g = 1 GeV. 

Figure IH also shows R and T for the same set of parameters as for Fig. [21 except 
that the magnetic field in the symmetric phase is B — 100 GeV 2 . 

Notice that in the cases depicted in Figs. [21 and El the reflection and transmission 
coefficients add up to unity, which in the present description of the scattering problem 
corresponds to the conservation of incoming flux of particles. 

To understand such flux conservation, let us recall that, according to Eq. Q, an 
incoming positive-energy particle from the symmetric phase encounters an effective 
potential barrier at the interphase of magnitude 

A V ~ ‘le'^B’(l + s' + 1/2) + /i 2 — ‘ItQB(l + s + 1/2), (34) 

which, for given quantum numbers of the incident flux, l and s, increases with the 
radial quantum number s' of the transmitted wave. Indeed, a one-dimensional version 
of the problem shows that such flux conservation, notwithstanding the relativistic 
context, is related to a potential difference that is mass-like. As argued below, this 
simulates the effect of a non-relativistic Schrodinger equation. We also compare this 
situation to the case of a gauge-like potential in the time component. 
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Figure 2: Reflection and transmission coefficients for an incoming wave function with 
quantum numbers s — 0, l — 1. The magnitude of the hypermagnetic held in the 
symmetric phase is B — 0.01 GeV 2 . The magnitude of the magnetic held in the 
broken-symmetry phase is B' = cos 9\yB. The coupling constant in the symmetric 
phase is e 0 — g' — 0.344, whereas in the broken-symmetry phase it is e' 0 = e = 
g' cosOw- The mass in the broken-symmetry phase is fi — 1 GeV. 

6 Toy models 

The conservation of hux can be understood with a simple one-dimensional model 
where the wave function satishes the Klcin-Gordon equation in two regions with 
different potentials. With the same conventions for regions / and II as in Sec. H3 we 
consider 


[d 2 - d‘l + Vf]vb = 0, z<0 

[d 2 t - d\ + V 2 2 + /i 2 ] T = 0, ^>0, (35) 

where Vf, V -2 are constant potentials for regions / and //, respectively. For E >V\ and 
E > \JV 2 + p 2 we have propagating incident, rehected, and transmitted momenta 
ki = \JE 2 — y 2 , k r = —kij q t = \JE 2 — V 2 2 — /i 2 , respectively. Using the explicit 

form for the conserved current given in Eq. H2SD, and assuming an incoming hux of 

positive-energy particles from region /, the reflection and transmission coefficients are 
given by 


{k,j - q t )(kj - q t )* 
{ki + q t ){ki + qt)* 1 
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Figure 3: Reflection and transmission coefficients for an incoming wave function with 
quantum numbers s — 0, l — 1. The magnitude of the hypermagnetic held in the 
symmetric phase is B = 100 GeV 2 . The magnitude of the magnetic held in the broken- 
symmetry phase is B' = cos 9\yB. The coupling constant in the symmetric phase is 
e 0 = g' = 0.344, whereas in the broken-symmetry phase it is e' 0 — e — g' cos 9w- The 
mass in the broken-symmetry phase is fi — 1 GeV. 


T 


4fc, 2 

(hi + q t ){ki + q t y ’ 


(36) 


these correspond to the ratio of the rehected (r) and transmitted (t) currents to the 
incident one, along the 5 axis 


jzr / jzi R 

izt/jzi = T R e(q t )/ki, (37) 

where the normalization is taken to the same density jo as in Eqs. (BED for all cases. 
As shown in Fig. 01 we get current conservation in this case, which we associate with 
the Schrodinger form of Eq. (j23)- This is not the case for a gauge-like potential in 
the time derivative, as in 


[(id t - Ri) 2 -d 2 ]T = 0, z < 0 

[(id t - V 2 ) 2 - d\ + p 2 ] T = 0, z > 0, (38) 

where Vj, V 2 are constant potentials for regions / and II, respectively. For E > Vj 
and E > V 2 + fi we have propagating incident, reflected, and transmitted momenta 
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Figure 4: Reflection and transmission coefficients in the one-dimensional with mass¬ 
like one-dimensional potential in Eq. f)38[) . with V± = 2.6 GeV in region z < 0, and 
V 2 = 2.3 GeV in region z > 0. The reflection and transmission coefficients add to 
one. The mass in the z > 0 region is fi = 1 GeV. 

hi = E — V 1 , k r = —hi, qt = y/(E — V 2 ) 2 + /i 2 , respectively. The currents are 

jzr/jzi (fH 

jzt/jzi = % - 17 T R e(q t )/ki, (39) 

Joj — V 2 

where T, R are still those of Eqs. m- It follows that within the positive-energy 
solution approximation used, the flux conservation ;j xl = j xr + j xt is valid only in 
the non-relativistic limit, and the ultra-relativistic limit E 00 , as shown in Fig. 0 
The flux non-conservation in the intermediate regime is associated to particle-pair 
creation. The trouble comes from the relativistic normalization in Eq. (HP . I 11 

practice, this means that a particle encounters a potential that makes possible the 
opening of additional channels contributing to the scattering processes, such as pair 
production. 

Such anomalous behavior is already known to arise in the context of the Dirac 
equation in the treatment of scattering and is called Klein paradox m- This signals 
the need to account for negative-energy solutions. In fact, the gauge-like potential in 
Eqs. m an d the current in Eq. (E2H1) produce the different current terms in Eq. (EH- 
It follows that the hypermagnetic-field problem resembles the Erst, Schrodinger-like, 
case. 
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Figure 5: Reflection and transmission coefficients with gauge-like one-dimensional 
potential in Eq. (1351) . and V\ = 2.6 GeV in region z < 0, and V 2 = 2.3 GeV in 
region z > 0. The reflection and transmission coefficients add to one only in the 
non-relativistic and ultra-relativistic limits. The mass in the z > 0 region is g = 1 
GeV. 


7 Conclusions 

In this work, we presented the exact solution for the scattering problem of scalar 
particles propagating within two regions of different but constant magnetic fields, 
separated by a planar interphase. The fields are assumed to be perpendicular to 
the interphase, while the solution for arbitrary directions is briefly explained in the 
Appendix. The two regions are defined by the condition that the particle is massless 
on one side of the interphase and massive on the other, such as during a first-order 
EWPT, outside and inside the nucleated bubbles, respectively. 

I 11 the general case in which the direction of the magnetic field is not perpendicu¬ 
lar to the wall interphase, we can decompose it into normal and parallel components. 
Recall that the normal component of the external field is continuous across the in¬ 
terphase and also that charged-particle trajectories consist of an overall displacement 
along the field lines superimposed to circular motion around the field lines. The 
relevant physical effect of the magnetic-field component directed along the wall is 
its ability to capture low momentum particles upon scattering since these can make 
transitions to states that can remain orbiting hied lines in this direction close to the 
wall. 

I 11 order to estimate the magnitude of such fields able to capture particles on the 
wall, let us look at the following classical argument: Equilibrium between the Lorentz 
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force and the centrifugal force gives for the radius of the orbit for a particle trapped 
in the wall R = p/(eB) where p is the particle’s momentum, e is its charge and B 
the strength of the magnetic held. Considering the physical situation of a wall with 
a finite thickness A, let us take R ~ A. Also, let us parameterize the strength of the 
magnetic held as B = bT 2 , where b is a dimensionless factor and T is the temperature, 
then b ~ p/(e\T 2 ). Taking the momentum of a typical particle to be of order of the 
temperature T, namely p ~ T and since A ~ T^ 1 (see for example Ref. 0), we 
obtain b ~ 0(10) which is already a very high value of the magnetic held. Although 
uncertainties on the strength of the magnetic helds at the EWPT are large, it is safe 
to assume that 6 10 (see for example Ref. |5], and references therein). 

We see that the loss of hux in these trapped modes for small held strengths is 
small, and therefore that the relevant magnetic-field component for the computation 
of transmission and reflection coefficients is the normal one, thus the magnetic held 
chosen in this paper. It is also worth mentioning that for a wall with non-zero thick¬ 
ness, realistic boundary conditions for the external gauge held imply the existence of 
a magnetic-field component along the direction of the wall and localized within the 
dimensions of the wall (see the discussion in Sec. IV of the second of Ref. [7]). In 
this case, similar arguments apply and the conclusion is that only very intense helds 
are able to produce an effect that involves a component of the magnetic held parallel 
to the wall. 

The problem at hand is solved through continuity conditions on a surface, rather 
than a point. The solution involves the expression of the wave function on one side 
of the interphase as a superposition of the complete set of solutions on the other side, 
and the continuity of the wave function and its derivative across the interphase. In 
order to numerically solve the problem, we resorted to a truncation of the system of 
equations involving the coefficients of the expansion. 

We found that in such a scheme, the flux of particles is conserved, which we 
associate with the Schrodinger nature of the magnetic-field equation. We showed in 
a simple one-dimensional model that such flux conservation can be understood as 
arising from the contribution of an effective mass-like potential barrier associated to 
momentum and radial modes, while, in the case of a gauge-like potential in the time 
component, the Klein paradox m, signals the need to account for negative-energy 
solutions. 

As previously found the scattering of chiral fermions during a first-order 

EWPT gives rise to an axial asymmetry that can be subsequently converted into a 
baryon asymmetry. The description of the scattering process used simplifying as¬ 
sumptions for the reflected and transmitted states. It will thus be interesting to 
explore the consequences of the exact method introduced in this work in the context 
of the Dirac equation on the generated axial asymmetry. This is work under progress, 
and will be reported as a sequel of the present one. 
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Appendix: Solution for arbitrary magnetic fields 

The method described in the paper can applied for the general case of two re¬ 
gions separated by a wall with fields Bi, Bo that are neither necessarily parallel to 
each other, nor perpendicular to the wall. Now we sketch the solution to the same 
scattering problem as above. 

We put the coordinates (x, y, z ) so that the wall separating the two regions is on 
the xz plane, so that y < 0 for region /, and y > 0 for region II (unlike the above 
case.) The fields define the plane B X B 2 perpendicular to the Bi x B 2 direction. We 
also choose the z axis along the line of the xz (wall) and B\B 2 planes intersection. 
The inclination of plane B\B 2 with respect to plane xz is measured by an angle <p 
drawn by rotation over the z axis. The directions of the B 1; B 2 fields on the BiB 2 
plane (which by definition lie on it) are measured by angles 6 and 9' between the z 
axis and the B 1; and B 2 directions, respectively. Indeed, the system is such that the 
spherical coordinates (B x ,9,(p) and {B 2 ,6',<p) also define the fields. Actually, each 
held defines another system of cylindrical coordinates such that travelling waves can 
be simply described. These are (r\,<pi,zi) and (r 2 ,cp 2 ,z 2 ), where the z x , z 2 axes are 
parallel to Bi, B 2 , respectively. 

Assuming an incoming wave /^ s (r 1; (pi, z x ) as in Eq. (TH2I) . it will be reflected to 

01, zi) = n.0,,z,) (40) 

nl 

(note that for didactic purposes, the arrangement is different than in Eq. (111(1 .1 and 
transmitted to 


*i r (r 2 ,02, z 2 ) = £ C„,,,/"(r 2 ,02, 22 ), (41) 

n'V 

respectively, on each side of the wall. Here we also sum over the angular momentum 
quantum number, which is no longer conserved, while the energy is. The conti¬ 
nuity of the wave function on the xz plane requires calculation of the overlap of 
the wave-function components. This is carried out transforming the coordinates 
{ri{x,y,z),(pi{x,y,z),z x {x,y,z)) and (r 2 {x,y, z),<p 2 (x,y, z), z 2 (x,y, z)) that are ob¬ 
tained respectively by the rotations. The condition for the wave function on the wall 
is y = 0, so that the overlapping amplitudes as in Eq. ITSlcan be calculated. We End 
that evanescent waves also play an important role as they provide the complete basis 
that allows for boundary conditions to be satisfied on a surface. 
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As for the case treated in the paper, one can use this method for an arbitrary 
incoming function, expanding it in the above basis. 
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